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Abstract 

We consider an exact Foldy-Wouthuysen transformation for the Dirac spinor field on the com- 
bined background of a gravitational wave and constant uniform magnetic field. By taking the 
classical limit of the spinor field Hamiltonian we arrive at the equations of motion for the non- 
relativistic spinning particle. Two different kinds of the gravitational fields are considered and in 
both cases the effect of the gravitational wave on the spinor field and on the corresponding spinning 
particle may be enforced by the sufficiently strong magnetic field. This result can be relevant for 
the astrophysical applications and, in principle, useful for creating the gravitational wave detectors 
based on atomic physics and precise interferometry. 
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I. INTRODUCTION 



The study of the gravitational effects on quantum-mechanical systems represents an im- 
portant issue, since all the physical objects, on both macroscopic and microscopic scales, 
are indeed located in a curved spacetime or in non-inertial reference frames. This sub- 
ject is currently of a special interest (see, e.g., an overview [1]). Despite the weakness of 
the gravitational interaction, its effects were actually observed at the quantum-mechanical 
level. In particular, one can mention the famous Colella-Overhauser- Werner (COW) ^ 
and Bonse-Wroblewski [3] experiments, in which the quantum-mechanical phase shift due 
to the gravitational and inertial forces was measured, thereby verifying the validity of the 
equivalence principle for the non-relativistic neutron waves. 

Among the different configurations of thegravitational field, the case of the gravitational 
waves appears to be especially interesting [4J]. This is related to the fact that considerable 
efforts are applied to the experimental search of the gravitational waves (see, e.g., Q]). 
Therefore, it may be important to achieve better understanding of the behavior of the 
quantum fields and particles under the action of the gravitational wave. 

In this paper, we study the dynamics of the Dirac particle in a plane gravitational wave. 
The plane-fronted gravitational waves represent an important class of exact solutions which 
generalize the basic properties of electromagnetic waves in flat spacetime to the case of 
curved spacetime geometry. The relevant investigation of the gravitational waves in general 
relativity has a long and rich history, see, e.g., [^j]. On the other hand, the approximate 
plane gravitational wave solution, that arises from the linearized Einstein's field equations, 
is usually used for the analysis of the various physical aspects related to the radiation, 
propagation, and detection of the gravitational waves. We will consider both cases, the 
approximate weak wave and the exact nonlinear plane-fronted wave. 

The general framework for the Dirac theory in curved spacetime was developed in many 
publications, see the early reference 0] , and the overview and the reference given in a recent 
paper We will use the notation and conventions of the latter work. 

The paper is organized as follows. In section [III we construct the universal form of an exact 
Foldy-Wouthuysen transformation for Dirac spinors and apply this scheme to the case of the 
background approximate plane gravitational wave and electromagnetic field. Several simpler 
particular (previously well known) cases are considered in Appendix. In section III 1 1 we use 
the result of the previous section to consider the particle Hamiltonian and corresponding 
equations of motion for the non-relativistic spinning particle interacting to the combined 
background of gravitational wave and electromagnetic field. In section IIVI we provide the 
comparison with the exact gravitational wave case and in the last section [V] we draw our 
conclusions. 
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II. LINEAR GRAVITATIONAL PERTURBATIONS AND UNIVERSAL FORM 
OF AN EXACT FOLDY-WOUTHUYSEN TRANSFORMATION 



Let us start from the relatively simple case of the usual (linear) gravitational waves. 
Using this case as an example, we shall also develop the general form of the exact Foldy- 
Wouthuysen transformation which can be applied to many previously explored cases and 
also to the more complicated case of the exact plane-fronted gravitational wave. 

The metric of the weak gravitational wave reads (see, e.g. [9[]) 

9ij = Vij + h ij , (!) 

where Tfcj = diag(— 1, 1, 1, 1) is the flat Minkowski metric and the nonzero components of 
the gravitational perturbation h^ v are (in the Cartesian local coordinates x % = (t,x,y,z)) 

^yy h>zz 2?7 , hyz h>zy . (2) 

Here v = v(ct — x) and u = u(ct—x) are the two functions which describe a wave propagating 
along the x axis. It is assumed that the functions v = v(ct — x) and u = u(ct — x) are small 
such that the linear approximation in these functions is valid. As usual, the gravitational 
wave can have 2 polarization states, and each of the functions v and u correspond to one of 
the possible polarizations. 

For the formulation of the Dirac theory in curved spacetime, we need the tetrad fields. 
The coframe 1-form reads: 

= c dt, tf 1 = dx, (3) 
$ 2 = (1 + v )dy + udz, (4) 
$ 3 = (1 - v )dz + udy. (5) 

The corresponding inverse vector frame is described by 

e = ^d t , e x = d x , (6) 
e 2 = (1 - v)d y - ud z , (7) 
e 3 = (1 + v)d z - ud y . (8) 

It is straightforward to construct the Riemannian connection (the Christoffel symbols). The 
nonzero components of the local connection 1-forms read: 

r 2 ° = T 2 = v'd 2 + u'd\ (9) 

r 3 ° = r 3 = u>$ 2 -v'$ 3 , (io) 
ry = -t 2 = v 'd 2 + u'd\ (ii) 
iy = -iV = u'$ 2 -v'$\ (12) 
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The primes denote the derivatives w.r.t. the argument of the functions: v' = dv(z)/dz,u' = 
du{z)/dz. 

One can check that the Cartan structure equation is fulfilled. Indeed, we have dd a + 
r /3 a A d 13 = 0, which demonstrates that this is the torsion-free Riemannian connection. 

Let us construct the Dirac operator on the background of the metric described above. 
The spinor covariant derivative is 



D a = e a \D, D:=d+-a af3 Y aP . 



We have a a/3 = 27^7^, and consequently, 



O P l a/3 



\ (7V r 20 + 7 3 7° r 30 + 7 V r i2 + 7V r 3 i) 



We shall use the standard lOj representation for the gamma matrices: 

I \ / a a 



7° = /? 



-I 



7 



-o a 



, a = 1,2, 3, 



where the Pauli matrices are: 



a 



1 

1 



a' 



-i 

1 



a" 



1 
-1 



The alpha matrices and the spin matrix are defined, as usual, by the relations 

if = 



a = 7 7 



-& 

CT> 



w 

o> 



Using these definitions, and substituting (!9|)- (fT2|) . we find 



a/3 



T a/3 = X - (c? + i E 3 ) (V 2 + u' 3 ) +^(a 3 -z E 2 ) (V tf 2 - v' & 
The Dirac operator is constructed as 7^-D^, and correspondingly, we obtain 



7%J f i ^ r Q , 



v' 
~2 



u' 
+ 2 



( 7 V-7 3 « 3 )+z( 7 2 S 3 + 7 3 S 2 
(7V + 7V) + i ( 7 3 E 3 - 7 2 E 



A direct computation shows that 

2 2 33 23,32 2v-*3 1 3v-*2 3^3 2\->2 n 

7 a — fa = 7 a + 7 a = 7 E + 7 E = 7 E — 7 E = 0. 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



Hence, all the terms on the r.h.s. of the equation ( fl9i) do vanish . As a result, the contribution 
of the spinor connection drops out completely and the Dirac operator is reduced just to the 
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partial derivative terms 



7 M Z^ = 7%Jd = - 1 °d t + 1 1 d x +\(l-v) 1 2 -ur 



+ 



(1 + t;) 7 3 -«7 : 



a 



Pl-d t + a l d x + [(1- 



a 2 — « a 3 



a, 



+ 



'1 + v ) a 3 — u a 2 



d z ). 



Now we are in the position to consider the covariant Dirac equation 

(ih^D^ - mc) ip = . 
After we recast it into the familiar Schrodinger form, this equation reduces to 

with the Hamilton operator of the form 

7i = mc 2 /3 + ca 1 p 1 + (1 — v ) ca 2 — u ca 3 p 2 + (1 + v ) ca 3 — u ca" 



p 3 , 



(21) 



(22) 



(23) 



(24) 



where p 1 , p 2 , p 3 are components of the momentum vector . 

The crucial observation is that this form of the Hamiltonian 
which admit the anticommuting involution operator, cf. 

BEL 

J = ry 5 /3. 



al 
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s into the class of models 



(25) 



The latter is Hermitian, = J, and unitary, JJ^ = J 2 = 1, and it anticommutes both with 
the Hamiltonian and with the (3 matrix: 



JH + HJ = 0, Jf3 + f3J = 0. 



(26) 



Consequently, the exact FW transformation can be constructed in this case. 

Moreover, if one introduces the (in the simplest case, constant and uniform) magnetic 
field, this important feature is not destroyed, since the only thing which is technically needed 
is to replace the operators of momentum by the expressions J? — - A . After performing 
this operation we arrive at the following Hamiltonian: 



H = f3mc 2 - ihc\a l (d x + —A x ) 
L inc 



(27) 



6 6 

+ (a 2 - va 2 - ua 3 )(d v + —A v ) + (a 3 - ua 2 + va 3 )(d z + —A z 

inc inc 



For a while we do not impose restrictions for the vector potential. Let us rewrite the 
expression (|27|) using new notations, which prove useful in what follows 13| 



H = (3mc 2 + a b K£d a + a a g a 



(28) 
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where 



/I 



Kl 



—ihc 







1 — v — ua 2 a 3 
VO 






1 + v — ua 3 a 2 J 



(29) 



9a 



-e {A l , A 2 - vA 2 - uA 3 , A 3 - uA 2 + vA 3 



(30) 



According to the standard prescription (the structure of exact FW for the scalar fields 
with non- minimal coupling to gravity has been discussed recently in 14j ) the fist step in de- 
riving the exact FW is to calculate the square of the Hamiltonian H 2 . In order to accomplish 
this, we define, additionally, conjugated quantities 



Ki = -ihc 



(1 



Vo 



o 

v + ua 2 o? 








(31) 



1 + v + ua at 



3^,2 



9a = {ot l gi , a 2 g 2 , a 3 g 3 ) , (32) 

and 

e = [ai€ , , «3e ) . {66) 

Using these notations, a direct calculation gives, after some algebra, the following square 
of the Hamiltonian H 2 : 

H 2 = m 2 c 4 + K^K h c d a d h + Kig b d a + g h K a b d a + g 2 

+ K^(d a K b c )d b + iY, d e bcd Kf{d a K:)d e + iZ^KfaPa 

- ie bcd Z d K a b g c d a + iZa^K^d^ + Kf(d a 9 b ) ■ (34) 

Here the partial derivatives inside the parenthesis act only on the content of these parenthe- 
sis. Let us notice that the components of the matrix K b depend exclusively on the metric 
while the components of the vector g a depend on both the metric and the electromagnetic 
potential A a . The expression (1341) is rather general and may be used not only for the linear 
gravitational waves but also in various other particular cases. In order to illustrate this fact, 
we consider several known cases in the Appendix. 



III. PARTICLE HAMILTONIAN AND EQUATIONS OF MOTION 

Starting from the equation (1341) . for the sake of simplicity, we consider only one polariza- 
tion of the gravitational wave. Namely, we choose u = 0. Then, we obtain 

K a h =Tq = -ihc(5 a b +T b a v) , g b = -e{52 + T b a v)A a , (35) 
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where 

T b a = I -1 I , (36) 

Vo o 1/ 

and there is no need anymore to use the notations 'g b and e abc . Thus we arrive at the 
Hamiltonian 

1 




H ~ 2^ P ^ + 2Tabv) ^ ~ ~ CAb) ~ ehc£ cadZ d d C (A b ) 

h 



+ -^(3e abc d a (v)Z c T b d (cp d -eA d )+{3mc 2 . (37) 
2mc 

A. Nonrelativistic limit 

The next step is to present the Dirac fermion ip in the form 

i) = ( £ j e ~ im(?t ' n , (38) 

and use the equation 

ihdtip = Hip (39) 

to derive the Hamiltonian for the 2-spinor ip. Inserting ([581 into (|39|) . we obtain the two- 
component equation 





-mc 2 + H) • (40) 



Using the fact that Hamiltonian is an even function, we obtain, in the (p sector, the Hamil- 
tonian 

H = - l —(5 ab + 2T ab v) \{cp a - eA a ){cp b - eA b ) - ehce cad a d cfA l 
2mc z 1 

+ ^- c e abc d a (v)a c T b d (cp d - eA d ) . (41) 

The r.h.s. of the last equation (Pill) is proportional to 1/m, that is the same level of the 
nonrelativistic approximation which one meets in the case of Pauli equation. Therefore one 
can expect that the same equation can be obtained starting from the original eq. ( 1281) . This 
calculation represents an efficient check of our results, so it is worthwhile to perform it now. 



Let us apply the standard procedure for deriving the Pauli equation (see, e.g. 15]). Using 
the representation (I4TJT) in (1281) and (139|) we apply the low-energy regime , that is assume 
that the term mc 2 is dominant (|mc 2 x| ^> |«7ic^x|). In this way we arrive at the equation 

lh Wt^ = 2^ {abKbda + ^SaK^IQdc + o°g B )<p . (42) 
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After some transformations, the r.h.s. of the last equation coincides with the r.h.s. of our 
eq. flUD. 

It is easy to check that, if the gravitational wave is absent v = 0, one obtains 
d 1 -> -> 

in dt p = ^w [1?(c1? ~ eA) ' ^ (c1? ~ eA ^ if ' (43) 

and finally, after a some algebra, a usual Pauli equation 

ih ® = [J-ff - e --tf - . -fifo . (44) 

dV [ 2m yF c ' 2mc ]y y 1 

B. Nonrelativistic spinor particle on the background of the gravitational wave and 
electromagnetic field 

One can use the the result ( l37j) for deriving the equation for the nonrelativistic spinor 



particle. Let us follow 16j, ll7j], where the similar calculation has been performed for the 



nonrelativistic equation on the background of torsion field. The same problem has been 



treated in [18( starting from the perturbative FW transformation. 

The Hamiltonian operador H corresponding to the energy ( 1371) is constructed in terms 
of the operators x a , p a and a a . The equations of motion have the form 

ih ~^f = ft*, h ) > ifc-Jjr = [P«» H \ ' ih ~^f = H \ ■ ( 45 ) 

In order to achieve the nonrelativistic limit one has to calculate the commutators of the 
operators in (145j) and take the limit h — » 0. Disrearding the 0(h) terms, we do not need 
to care about the ordering of the operators, and then finally we arrive at the semiclassical 
equations of motion for the spinning particle on the background of the gravitational field 
and constant magnetic field. The result can be presented using the notation A' a = T a bA b . 
Then the equations become 



dx a 

dt m 



-(5 ab + 2T ab v) (p b --A b ) , (46) 

71 \ C J 

T bc d a v ( Pb --A b ) ( Pc --A c ) + (S bc + 2T bc v) (p b --A b ) —d a A c ,(A7) 



da a e c 

—r~ — ZabcC 

dt mc 



^ + 2v rot (T)l " - — (cp d - eA d ) a e [T ed (d a v) - T ad (d e v)} . (48) 



mc 



Let us notice that the first equation here demonstrates deviation from the usual relation 
between the Lagrangian velocity and the canonically conjugated momenta. The difference is 
due to the gravitational wave. Taken together, the first two equations define the coordinate 
dependence of the particle. Using these two equations, one can construct the analog of 
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Lorentz force in the presence of gravitational wave 



m x a 



e 
c 



e ,„ _ , dA b 



{5 ab + 2T ab v) lit x ~£] b - -{5 ab + 2T ab v 



c at 



+ 2mT afe — x b - mT bc x b x c (d a v) . (49) 



The last equation in (14811 in describing the spin dynamics of the particle. It is remarkable 
that this equation depends on the velocity. Let us notice that the same property holds for 
ihe spinning particle in the external torsion field for both nonrelativistic 16j and relativistic 
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cases. 

Perhaps the spin dynamics of the particle is the most interesting result here. In order to 
understand this point, let us make the following observation. It is easy to notice that the 
last equation in (j4*8l) describes spin precession even for the case when the magnetic field is 
absent. The effect is due to the presence of the matrix T fc a defined in eq. (l36l) and the rot A'. 
The last vector is distinct from the magnetic field and maybe nonzero when B = rot A = 0. 
Therefore, there is a possibility to observe the spin precession without magnetic field and 
this can, in principle, become a new basis for the gravitational waves detector of the new 
type. We expect to consider this issue in more details in the near future. A detailed analysis 
of the definitions of the momentum and spin dynamical operators (along the lines recently 
done for the static gravitational field [20|) will be required. 



C. Vanishing electromagnetic field 

In the absence of the electromagnetic potential, A a = 0, the equations of motion can be 
intergrated. With x a = (x,y, z), the second and the third equations in (T49"]) read 

U = -2y%, (50) 
dt y y dt' y ' 

d . n . dv 1 . 

— z = 2z—. 51 
dt dt 

These equations are easily integrated, yielding the components of the velocities as the func- 
tions of v : 

y = Ye~ 2v , z = Ze 2v , (52) 

with the integration constants Y and Z. The integration of the equation for x is slightly 
more nontrivial. This equation reads 

x = (y 2 - z 2 ) d x v = f(v) d x v, (53) 

with f(v) = Y 2 e~ Av — Z 2 e 4v . Recalling that v = v(£) with £ = ct — x, we notice that 
d x v = —dv/dl;, and furthermore £ = —x. Accordingly, we can recast (153!) into the equation 
for £: 

2*=-|r. U= 1 -(Y 2 e~^ + Z 2 e^). (54) 
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Multiplying with £, we find the first integral 

+ u(0 = i, 

and the solution £ = is then obtained in quadratures 



U(0 



t-tn 



(55) 



(56) 



The explicit form of the solution £ = £(£) of course depends on the integration constant I and 
on the explicit form of the wave function v = i>(£) that will determine the "potential" £/(£). 
For the harmonic wave, for example, v = Vq cos£, and then Z7(£) = y 2 e~ 4cos £ + Z 2 e 4cos ^ 
(with some new constants F and Z). After finding we can use it for the final integration 
of the first order equations (I5"2j) . Thus, we finally obtain the coordinates of the particle as 
functions of time y(t), z(t) and x(t) — ct — 

In order to complete the analysis of the particle dynamics, we have to solve the equation 
for the spin. In the absence of the electromagnetic field, it reads: 

\Jv% + Vl (d x v) 



ni 



M ab a b , 



with the matrix 



( ix y -n z ^ 




-TTy 

7T Z 



Here 7r„ 



Pyl \jp 2 y + Pli n z = p z / y^| + pf. Note that p y = p z 



0, see 



(57) 



(58) 



The integration 



is then straightforward, yielding the final result for the dynamics of the spin 
!? = exp [s(i)M] • "^o , where s(t) = 



Pl + P 2 z 



rn 



dtdJv) 



(59) 



to 



Note that since for the cubic term we have M 3 = — M, the matrix exponential actually 
contains only the terms M and M 2 , and it can be written explicitly as 



a s(t)M 



1 + M sin s + M 2 (1 - cos s) 



The solutions fl52l) . fl56|) . fl59|) show that the nonrelativistic spinning particle has very 
peculiar behaviour in the field of the weak gravitational wave. 

IV. COMPARISON WITH EXACT GRAVITATIONAL WAVE 



The exact plane-fronted gravitational wave, in the simplest case is described in the Carte- 
sian local coordinates x l = (t, x, y, z) by the line element ds 2 = g i jdx l dx : ' with the metric 
(see 0,0,3) 



9ij = Vij + h. 



i.i ■ 



(60) 
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where the nonzero components are 

h tt = h xx = -U, h tx = h xt = U (61) 

expressed in terms of a function U(£,y,z). It can depend arbitrarily on £ = ct — x, and is 
a harmonic function in the two last variables, i.e. A( 2 )U = (dy y + d^ z )U = 0. The coframe 
reads 

= (l + ^ cdt - ^ dx, (62) 

tf 1 = ^ crft + (i - ^) (63) 

tf 2 = A/, tf 3 = d^. (64) 



The inverse frame is easily found: 



eo= (l-^)Ift-^, (65) 

ei = + ( 66 ) 

e 2 = d w , e 3 = d 2 . (67) 



From this we can verify straightforwardly that 



% 
4 

Accordingly, we find that 



\ I> = ( 7 ° - 7 1 ) \ ^d y U + ^d z U) {f - . (68) 



7%J (1 ^ r Q/3 ) = - (7° - 7 1 ) 2 ^ ( W + 7^) = 0, (69) 

since (7 — 7 1 ) 2 = 0. Thus, just like in the case of an approximate wave, the spinor connec- 
tion term drops out completely from the Dirac equation. 
Consequently, the Dirac operator has the form 

l a D a = ^e a \d = - c [l° + |(V - 7°)] ft + [7 1 + - 7°)] d x + ! 2 d y + ^d z . (70) 

It proves very useful, before we go to the exact FW transformation, we make a Lorentz 
transformation of the coframe — > = A a /3 ^ /3 , using the matrix (written in the 2x2 
block form with and 1 as the 2x2 zero and unit matrices, respectively) 

Under this transformation, the tetrad frame is changed to e a — > e' a = (A^ 1 ) /3 a e / 3, whereas 
the local Lorentz connection transforms to T'J 3 = A /3 u T fl L '(A' 1 )^ a + A /3 7 rf(A" 1 ) 7 Q ,. One can 
then verify that the transformed spinor connection term do contribute to the Dirac operator 

1%} ( J ** I"*) = (7° - V ) + jjj^j 7° (,Y <V + tVW) . (72) 
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where U' = dU/d^. Another contribution comes from the ordinary derivative term 



l a e' a \d = l Z 



U 



7 



9« + 7 2 9» + 7 3 ^- 



(73) 



Collecting all together, we can finally write the Dirac equation (122]) in the Schrodinger form 
63). As a last s^ep we perfan the rescaling transformation 0~the wave function 

= (1 - U)H . 

In this way we arrive at the final form of the Hamiltonian that is explicitly Hermitian: 



H' = f3mc 2 V + cp x - Ul - a 1 ) (V 2 p x + p x V' 



+ 



2 

ihc 



a 2 (V Py + p y V) + a 3 (Vp z + p z V) 



a 



a 2 d y {V) + a 3 d z (V) 



(74) 



Here V = l/y/l^U. 

Let us notice that switching the gravitational wave off with U = (hence V — 1), we 
recover the Hamiltonian of the free particle. The Hamiltonian (IT4"|) is not anticommuting with 
the matrix J = ij 5 [3. As a result, the condition for performing the exact FW transformation 
({J, H'} = 0) is not satisfied in this case. 



V. DISCUSSION AND CONCLUSION 



In this paper, we have derived an exact Foldy-Wouthuysen transformation for the Dirac 
spinor field on the combined background of a gravitational wave and constant uniform mag- 
netic field. The motivation for the presence of the magnetic field is to check the possibility 
of the amplification of the influence of the gravitational wave on a Dirac particle. According 
to our calculations (1521) . (|47|) and (|48j) such an effect is possible. This result can be relevant 
for the astrophysical applications and, in principle, could be useful for improving the gravi- 
tational wave detectors based on atomic physics and precise interferometry (see, e.g., recent 
discussion in (23)]). The actual dynamics of the Dirac particle in a plane gravitational wave 
will be studied separately. 
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APPENDIX A: FREE PARTICLE 



Let us consider several simple known particular cases when the general formula flMl) can 
be applied. 

For the case of a free particle, we have K% = —ihc8% and g a = 0, because u = v = 
and A\y = 0. Then 



H = /^m 2 c 4 + cV , (Al) 
that is a well-known known textbooks result. 

APPENDIX B: PARTICLE IN A MAGNETIC FIELD 

In this case K% = —ihc5% and g a = —eA a , because u — v — 0. The expression for the 
Hamiltonian is 



H = /3^m 2 c* + (cf -el) 2 -hcell.~B* . (Bl) 
This is exactly the result obtained by Eriksen and Kolsrud for this case. 

APPENDIX C: PARTICLE WITH ANOMALOUS MAGNETIC MOMENT IN A 
STATIC MAGNETIC FIELD 

In this case we start from the Hamiltonian with 

Kl = -ihc5 a b 

and g a = —eA a + a a ///E .B , and again with u — v — 0. 

This version is a bit more complicated because one has to account for the commutators of 
g a with a b and (3. After some calculations we arrive at the following form of the Hamitonian 

H 2 = m 2 c 4 + (cf - elf - 2(t I mc 2 ll.l$ + fxjB 2 + ^0 ■ (~B x - ~f x if) . (CI) 



Once again, this expression is in a perfect agreement with the result obtained by Eriksen 



and Kolsrud [11] 
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